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Abstract 

In a gas of ultracold atoms whose scattering length is controlled by a magnetic Feshbach reso¬ 
nance, atoms can be associated into universal dimers by an oscillating magnetic field. In addition 
to the harmonic resonance with frequency near that determined by the dimer binding energy, there 
is a subharmonic resonance with half that frequency. If the thermal gas contains dimers, they can 
be dissociated into unbound atoms by the oscillating magnetic field. We show that the transition 
rates for association and dissociation can be calculated by treating the oscillating magnetic field as 
a sinusoidal time-dependent perturbation proportional to the contact operator. Many-body effects 
are taken into account through transition matrix elements of the contact operator. We calculate 
both the harmonic and subharmonic transition rates analytically for association in a thermal gas 
of atoms and dissociation in a thermal gas of dimers. 
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I. INTRODUCTION 


The use of magnetic Feshbach resonances to control the interaction strengths of ultracold 
atoms has led to signihcant advances in our understanding of strong interactions in few-body 
and many-body physics. The effects of time-dependent strong interactions can be studied 
by using a time-dependent magnetic held. A particularly interesting case is a sinusoidally 
modulated magnetic held. Atoms can be associated into universal molecules composed of 
atoms with a large scattering length by modulating the magnetic held with a frequency near 
that determined by the binding energy of the molecule. The measurement of the binding 
energy of a molecule by the resonance in the oscillation frequency is called magnetic-field 
modulation spectroscopy or sometimes wiggle spectroscopy. 

Modulation of the magnetic held was pioneered by Thompson, Hodby, and Wieman to 
associate ®^Rb atoms into dimers [1]. Papp and Wieman used magnetic-held modulation 
spectroscopy to measure the small binding energies of dimers composed of ®^Rb and ®^Rb 
atoms [2]. Weber et ah used a resonantly modulated magnetic held to associate and 
®'^Rb atoms into dimers and to measure their binding energies [3]. They also observed 
subharmonic resonances at half the frequency determined by the binding energies of the 
dimers. Lange et ah used magnetic-held modulation spectroscopy to measure the binding 
energies of ^^^Cs dimers [4]. Pollack et ah used a modulated magnetic held to excite collective 
modes in a Bose-Einstein condensate of ^Li atoms [6]. Machtey et ah used a modulated 
magnetic held to associate ^Li atoms into Ehmov trimers [7]. Dyke, Pollack, and Hulet 
used magnetic-held modulation spectroscopy to measure the binding energies of ^Li dimers 
in both a Bose-Einstein condensate and a thermal gas [8]. In the thermal gas, they also 
observed a subharmonic resonance. Smith recently pointed out that a sinusoidally oscillating 
magnetic held near a Feshbach resonance can also be used to control the scattering length, 
and he showed that the resonance parameters are universal functions of the magnetic held 
[9]. 

A theoretical treatment of the association of atoms into dimers by an oscillating magnetic 
held was hrst presented by Hanna, Kohler, and Burnett in 2007 [10]. An alternative approach 
was recently developed by Brouard and Plata [11]. Both groups described the two-atom 
system by a two-channel model consisting of a continuum of atom-pair states and a discrete 
molecular state. They calculated the probability for the association of atom pairs into 
dimers as a function of time by solving the time-dependent Schrodinger equation for the two 
coupled channels. In Ref. [11], some qualitative aspects of the harmonic and subharmonic 
association processes were derived analytically. The results for association probabilities from 
both groups were completely numerical. 

A much simpler approach to this problem was recently introduced in Ref. [12]. It was 
inspired by Tan’s adiabatic relation, which expresses the change in the energy E of a system 
due to a change in the scattering length a in terms of an extensive thermodynamic variable 
that is conjugate to 1/a called the contact C [13]. In the case of fermions with mass m and 
two spin states, the adiabatic relation is 

d ^ 

d{l/a) Anm 

(In the case of identical bosons, the right side should be multiplied by 1/2.) Ref. [12] 
pointed out that the transition rates for the association of atoms into universal dimers 
can be calculated by treating the oscillating magnetic field as a sinusoidal time-dependent 
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perturbation proportional to the contact operator. The association rates were calculated for 
a thermal gas of atoms and for a dilute Bose-Einstein condensate of atoms. In this approach, 
many-body effects are taken into account through transition matrix elements of the contact 
operator. 

In this paper, we extend the approach of Ref. [12] to the dissociation rates of universal 
dimers and to subharmonic transitions. In Sections II, we derive general formulas for the 
harmonic and snbharmonic transition rates using time-dependent perturbation theory. In 
Sections III and IV, we calculate the leading harmonic contributions to the association 
rate in a thermal gas of atoms and the dissociation rate in a thermal gas of dimers. They 
come from a hrst-order perturbation in the contact operator. In Sections V and VI, we 
calculate the dominant snbharmonic contribntions to the association rate in a thermal gas 
of atoms and the dissociation rate in a thermal gas of dimers. They come from a second- 
order perturbation in the contact operator. In Section VII, we apply our results for the 
association rate to a thermal gas of ^Li atoms. In Section VIII, we summarize previous 
theoretical treatments of association into dimers using a modulated magnetic field, and we 
compare them with onr results for association. 


II. TRANSITION RATES 

In this section, we derive general formulas for transition rates at first order and second 
order in time-dependent pertnrbation theory. We focus on the case of fermionic atoms with 
eqnal mass m and two spin states. (We also give the corresponding results for identical 
bosons.) 

A. Perturbing Hamiltonians 

Near a magnetic Feshbach resonance, the scattering length a of the atoms is a fnnction 
of the magnetic held: 

a{B) = abg[l - A/{B - Bq)], (2) 

where Obg is the background scattering length and Bq and Bq + A are the positions of the pole 
and the zero of the scattering length, respectively. We consider a time-dependent magnetic 
held that has a constant value B for t < 0 and is modulated with a small amplitude b around 
the average value B for t > 0: 

B{t) = B t < 0, 

= B + bsm{ut) t > 0. (3) 

If the oscillating magnetic held is inserted into Eq. (2), it implies a time-dependent scattering 
length a(t). 

Tan’s adiabatic relation in Eq. (1) implies that the leading perturbation in the Hamilto¬ 
nian for t > 0 is proportional to the contact operator: 

ft 2 / 1 1 \ 

^ ^ ^ ^ 47rm \a{t) aj ' ^ ^ 

where a = a{B) is the scattering length in the absence of the modulated magnetic held. 
In Appendix A, quantnm held theory methods are used to argue that this is the only 
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perturbation that contributes in the zero-range limit. The inverse scattering length can be 
expanded in powers of b: 
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a{t) 
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(5) 


The coefficients of the powers of b have well-behaved limits as B approaches the Feshbach 
resonance at B^. Inserting the expansion in Eq. (5) into Eq. (4), we can identify terms in 
the perturbing Hamiltonian that are hrst and second order in b: 


HAt) = 


H,Jt) = 
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(In the case of identical bosons, the right sides shonld be multiplied by 1/2.) If |6| |A|, 

the effects of Hi and H 2 can be taken into account as time-dependent perturbations. The 
hrst-order perturbation in Hi drives transitions to states whose energies are higher or lower 
by huj. We refer to such transitions as harmonic transitions. The first-order pertnrbation in 
H 2 and the second-order perturbation in Hi both drive transitions to states whose energies 
differ by 0 or ±2^0;. We refer to transitions to states whose energies are higher or lower by 
2huj as suhharmonic transitions. 


B. Fermi’s Golden Rule 


We hrst consider transitions from the hrst-order perturbation in Hi. We take the initial 
state |f) to be an energy eigenstate with energy H*. We consider the transition to a distinct 
energy eigenstate |/) with energy Ef. At hrst order in pertnrbation theory, the probability 
amplitude for the hnal state |/) at time T is 



ih . 

( ) 

'^i{u>fi+u))T _ 


STTmaog 

{{A + Bo-By) 

OJfi -F OJ 

UJfi-UJ 


ifm), 


(7) 


where Ufi = {Ef — Ei)/h. The two terms inside the brackets have absolnte values that 
increase linearly with T in the limits ca/j —)■ —u and ojfi ^ + 0 ;, respectively. By applying 
Fermi’s Golden Rule, we obtain the transition rate summed over hnal states |/): 
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(In the case of identical bosons, the prefactor should be multiplied by 1/4.) This transition 
rate is non-zero only for hnal states whose energy dihers from Hj by ihca, so it contributes 
to the harmonic transition rate Fi(a;). The superscript (1) on F^^^ indicates that it comes 
from the hrst-order perturbation in Hi. 

We next consider transitions from the hrst-order perturbation in H 2 . At hrst order in 
perturbation theory, the probability amplitnde for the hnal state |/) at time T is 
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Inside the brackets, we have shown explicitly only those terms whose absolute values increase 
linearly with T in the limits ca/j —)■ ±20;. By applying Fermi’s Golden Rule, we obtain the 
transition rate summed over hnal states |/): 



2567r±n^ \ {A + Bq - By 


f ± 


( 10 ) 


(In the case of identical bosons, the prefactor should be multiplied by 1/4.) This transition 
rate is nonzero only for hnal states whose energy differs from Ei by A2hu, so it contributes 
to the subharmonic transition rate r 2 (o;). The superscript (2) on T^ indicates that it 
comes from the hrst-order perturbation in H 2 . The subharmonic transition rate in Eq. (10) 
is determined by the same transition matrix element {f\C\i) of the contact operator as the 
harmonic transition rate in Eq. (8). It can be expressed in terms of the leading harmonic 
transition rate at twice the frequency: 


p(2) 


(cn) 


1 / b 
4\A + Bo-B 


r«(2cn). 


( 11 ) 


Finally we consider transitions from the second-order perturbation in Hi. At second 
order in perturbation theory, the probability amplitude for the hnal state |/) at time T is 
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where the sum is over intermediate states \m) distinct from |i). Inside the brackets, we have 
shown explicitly only those terms whose absolute values increase linearly with T if ujfi is 
±2ci;. By applying Fermi’s Golden Rule, we obtain the transition rate summed over hnal 
states I/): 
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(13) 

(In the case of identical bosons, the prefactor should be multiplied by 1/16.) This transition 
rate is nonzero only for hnal states whose energy dihers from Ei by ±2ha;, so it contributes 
to the subharmonic transition rate F 2 (a;). The superscript (1,1) on indicates that it 
comes from the second-order perturbation in Hi. There is an additional factor of l/Obg in 

the prefactor for compared to F^^^. The relative importance of these two contributions 
is determined by the canceling length scales provided by the contact matrix elements and the 
frequency denominator. If and F^^^ have the same order of magnitude, the interference 
between the hrst-order perturbation in H 2 and the second-order perturbation in Hi would 
have to be taken into account. By explicit calculations of subharmonic transition rates in a 
thermal gas, we will hnd that the additional dimensionless factor in is (a/obg)^- Thus 
F 2 ^’^^ is much larger than F^^^ if B is near a Feshbach resonance. 
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C. Thermal System 


The transitions rates in Eqs. (8), (10), and (13) apply to an initial state |i) that is 
an energy eigenstate. A thermal system is described instead by a density matrix. For 
a completely thermalized system, the density matrix is p = exp(—/3if)/Tr(exp(—/3iJ)), 
where H is the Hamiltonian and [3 = l/ksT. By expressing the modnlus-sqnared of an 
amplitnde as the prodnct of the amplitnde and its complex conjngate, the dependence on 
the initial state in Eqs. (8), (10), and (13) can be pnt in the form of the projection operator 
mnltiplied by a fnnction F{Ei) of the initial energy that inclndes the freqnency delta 
fnnction. If the density matrix p is diagonal in an energy basis, the transition rate is obtained 
by making the snbstitntion 

r(Bj)|*)(i|—>5^f(B.)|*>(*|p|i){i|. (14) 


D. Homogeneous System 

The contact operator C is an extensive variable. It can be expressed as the integral over 
space of the contact density operator: 



If the initial and hnal states are homogeneous systems, we can simplify the transition rates 
by expressing them in terms of matrix elements of the contact density operator. 

The harmonic transition rate r^^^(ci;) in Eq. (8) and the subharmonic transition rate 
r 2 ^^(a;) in Eq. (10) involve the factor |(/|C'|i)p. By inserting the expression for C* in Eq. (15), 
we obtain matrix elements of the contact density at two different positions. We can use 
translational invariance to put both operators at the same position r. One of the integrals 
over space then gives a momentum-conserving delta function. The resulting expression for 
the modulus-squared of the transition matrix element is 

= |(/|CWI*)|t ( 16 ) 

/ / 

where Ki and Kf are the total wave vectors of the initial and hnal states of the homogeneous 

I 12 

system, respectively. Homogeneity implies that (/|C(r)|i) is independent of the position 
r. Thus the integral Jd^r in Eq. (16) just gives a factor of the volume V. 

The subharmonic transition rate r 2 ^’^^(ci;) in Eq. (13) involves the product of four matrix 
elements of C. By inserting the expression for C in Eq. (15), we obtain matrix elements 
of the contact density at four different positions. We can use translational invariance to 
put all four operators at the same position r. Three of the integrals over space then give 
momentum-conserving delta functions. The resulting expression for the factor in Eq. (13) 
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that involves matrix elements of C is 
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where Km and Km' are the total momenta of the intermediate states \m) and |m'), respec¬ 
tively. Homogeneity implies that the integrand is independent of the position r. Thus the 
integral J d^r just gives a factor of the volume V. 


E. Local Density Approximation 

For a many-body system whose number density varies slowly with the position r, the 
transition rate can be simplihed by using the local density approximation. The transition rate 
is an extensive quantity. For a homogeneous system, the expressions obtained by inserting 
Eq. (16) or Eq. (17) into the transition rate have an explicit factor of the volume J d^r = V. 
If the transition rate is also proportional to the total number iVj of some type of particle, 
the additional factor must be the intensive combination Ni/V. For a homogeneous system 
consisting of fermionic atoms with spin states 1 and 2, the association rate is proportional 
to N 1 N 2 . The local density approximation for the association rate in a system with local 
number densities ni{r) and n 2 {r) can be obtained by making the substitution 

N 1 N 2 /V —)■ J d^rni{r)n2{r). (18) 

For a homogeneous system consisting of dimers, the disssociation rate is proportional to their 
total number Ad. The local density approximation for the disssociation rate in a system 
with local number density nBir) can be obtained by making the substitution 



III. HARMONIC ASSOCIATION RATE 

A pair of atoms with a large positive scattering length can be associated into a universal 
dimer by an oscillating magnetic field. In this section, we calculate the harmonic association 
rate in a thermal gas of atoms. We also give the subharmonic association rate from first- 
order perturbation theory. We consider a gas of atoms that is in thermal equilibrium at 
temperature T. For simplicity, we take the number densities ni and n 2 of the atoms to be 
sufficiently low that their distributions are given by Boltzmann statistics instead of Fermi- 
Dirac statistics. 
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A. Initial and Final States 


We first consider a homogeneous gas consisting of Ni atoms of spin state 1 and N 2 
atoms of spin state 2 in a volume V. The two spin states interact with a large positive 
scattering length a. The universal dimer has a small binding energy h?'/ma?. For a gas in 
thermal equilibrium, the harmonic transition rate is given by Eq. (8) with the substitution 
in Eq. (14), where p = is the density matrix for the thermal gas of atoms. To simplify 
the presentation, we will temporarily ignore the frequency delta function, which depends on 
the energy Ei of the states in the density matrix. The terms in Eq. (8) that depend on the 
contact operator can then be expressed compactly as X]/(/I^PgasC'l/)- We will insert the 
frequency delta function at the end of the calculation. 

In the low-density limit where 3-body and higher-body correlations can be neglected, the 
density matrix p^^s can be expressed in terms of the density matrix ppair for a pair of atoms 
in thermal equilibrium: 

Y.(f\Cp„.C\!) = iV,iV2^{/|Cpp.„C|/>. (20) 

/ / 


The factor N 1 N 2 is the number of pairs of fermions in the two spin states. (For a gas of 
N identical bosons, the number of pairs is N‘^/2.) The pair density matrix is normalized: 
Tr(ppair) = 1- On the left side of Eq. (20), the sum over / is over many-body final states 
that include a single dimer. On the right side, the sum over / is over two-atom final states 
that consist of a single dimer. The density matrix for a pair of atoms in thermal equilibrium 
is 


Ppair = ^ j j exp(-/ 9 fiVl'^/ 4 m - fc)(fc,/f I, 


( 21 ) 


where (5 = l/ksT and At is the thermal deBroglie wavelength for an atom with mass m: 


At — 


\J 2nh‘^ /mk-QT. 


( 22 ) 


The two-atom states \K, k) in Eq. (21) are labeled by the center-of-mass wave vector K = 
ki + k 2 and the relative wave vector k = (fci — k2)/2. The integrals over the wave vectors 
are dehned by 



(23) 


The wave vector states have delta-function normalizations: {k'\k) = {27r)^S^{k' — k). In the 
case k' = k, the inhnite norm can be expressed as a factor of the volume: {k\k) = V. The 
energy of a pair of atoms in the state \K, k) is 


Eaa = /^m + h‘^k‘^/m. 


(24) 


The sum over hnal states on the right hand side of Eq. (20) can be expressed as an 
integral over the wave vector of a dimer: 


/ 


(/IC'ppairC'l/) = / (fcDlC'ppairC'lfco). 
J ku 


(25) 


The energy of the dimer is 


Ed = -\- fi^k^/Am. 


( 26 ) 



B. Matrix Elements 


Because the system is homogeneous, the analog of Eq. (16) can be used to express the 
contact operators C on the right side of Eq. (25) in terms of contact density operators at 
the same position r. The wave vector delta function in Eq. (16) reduces to S^{kD — K), and 
it can be used to integrate over fco- The frequency delta function in Eq. (8) reduces to 

27r(5((ED — Eaa)/^ ± ca) = 2ti5 (lo — h/ma^ — /m). (27) 

± 

In the sum over ica, only the +a; term contributes. 

The expression for the transition rate has been reduced to matrix elements of the con¬ 
tact density operator of the form {kY)\C{r)\K,k). The matrix element is calculated in 
Appendix B, and is given by Eq. (B12): 

{fcD|C(r)|/f,fc> = ^^^. ( 28 ) 

The Gaussian integral over K can be evaluated analytically. The sum over hnal states of 
the matrix element in Eq. (25) reduces to 

AT AT /“OO 

X](/|CpgasC^|/) = 128v^7raA^^^ y dk ^^ exp{-/3h‘^k‘^/m). (29) 

Before integrating over k, this must be multiplied by the frequency delta function in Eq. (27). 


C. Harmonic Association Rate 


Our hnal result for the harmonic association rate r[^^(ci;) in the homogeneous gas can 
be obtained from Eq. (8) by replacing |(/|G|i)p by the right side of Eq. (29), replacing 
the sum of frequency delta functions by the right side of Eq. (27), and then using the delta 
function to integrate over k. The local density approximation can be implemented by making 
the substitution for N 1 N 2 /V in Eq. (18). The threshold angular frequency for association 
is h/ma^: the emission of a smaller energy from a pair of atoms is not enough to allow a 
transition to dimer. For u > h/ma^, the harmonic association rate is 



(ca) 


where 


2V2h^ / bA 
w?al^a \{A + Bo- By 



Arp/T)(ta) ni-2 2 t \ I \ 

- - exp(—K^(a;)/m), 

ijJ 

(30) 


k(cj) 


a/ mu/h — \/o?. 


(31) 


(The harmonic association rate in a thermal gas of identical bosons with large scattering 
length was calculated in Ref. [12]. It can be obtained from Eq. (30) by replacing ni{r)n 2 {r) 
by n^(r)/2, where n{r) is the local number density of identical bosons.) If ksT /ma^, 
the harmonic association rate in Eq. (30) has a narrow peak with a maximum when u is above 
the threshold h/ma^ by approximately kBT/2h. For large frequency, the rate decreases as 
exp(—hcn/Zc^T). 
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D. First-Order Subharmonic Association Rate 


According to Eq. (11), the contribution r 2 ^^(ci;) to the subharmonic association rate from 
the hrst-order perturbation in H 2 can be expressed in terms of the harmonic association 
rate in Eq. (30) at twice the frequency. The threshold angular frequency for subharmonic 
association is \{h/ma^). For u > the subharmonic association rate is 



Am‘^al^a \{A + Bo-By) 



Arp/^(2cu) , 9/^ \ / \ 

—- exp(-/3h k^(2lo)/ m), 

CJ 

(32) 


where k(2(^j) is the function dehned in Eq. (31) with a; replaced by 2a;: 


k(2(^j) = ^/2mu/h — l/cB. (33) 

If ksT /ma?, this contribution to the subharmonic association rate has a narrow 

peak with a maximum when a; is above the threshold \{h/ma^) by approximately ksT/Ah. 
The height of the peak is smaller than that of the harmonic association rate by the factor 
[V(A + 5o-5)]V4. 


IV. HARMONIC DISSOCIATION RATE 

A universal dimer can be dissociated by an oscillating magnetic held into its constituent 
atoms. In this section, we calculate the harmonic dissociation rate in a thermal gas of 
dimers. We also give the subharmonic dissociation rate from hrst-order perturbation theory. 
We consider a gas of dimers in thermal equilibrium at temperature T. For simplicity, we 
take the number density ud of dimers to be sufficiently low that their distribution is given 
by Boltzmann statistics instead of Bose-Einstein statistics. 


A. Initial and Final States 

We hrst consider a homogeneous gas consisting of Nd dimers in a volume V. If the gas is 
in thermal equilibrium, the harmonic transition rate is given by Eq. (8) with the substitution 
in Eq. (14), where Pgas is the density matrix for the thermal gas of dimers. To simplify the 
presentation, we will temporarily ignore the frequency delta function, which depends on 
the energy Ei of the states in the density matrix. The terms in Eq. (8) that depend on 
the contact operator can be expressed compactly as X]/(/I^PgasC'|/)- We will insert the 
frequency delta function at the end of the calculation. 

In the low-density limit where correlations between dimers can be neglected, the density 
matrix pgas can be expressed in terms of the density matrix pdimer for a single dimer in 
thermal equilibrium: 


5I(/|CPs„C'|/) = No 5^(/|Cpdi„,C|/). (34) 

/ / 

The dimer density matrix is normalized: Tr(pdimer) = 1- On the left side of Eq. (34), the 
sum over / is over many-body hnal states that include an unbound pair of atoms. On the 
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right side, the sum over / is over two-atom hnal states that consist of an unbound pair of 
atoms. The density matrix for a dimer in thermal equilibrium is 


Pdimer 


Ay 



{-I3h‘^kl/Am) |fcD)(fcD|, 


(35) 


where At is the thermal deBroglie wavelength for an atom in Eq. (22). The energy Ed of 
the dimer is given in Eq. (26). 

The sum over hnal states on the right side of Eq. (34) can be expressed as integrals over 
the total wave vector and the relative wave vector of a pair of atoms: 

V(/|C'pdimerC'|/) = / / (K, fc|C'pdimerC'|K, fc) . (36) 

f J K Jk 


The energy Eaa of the pair of atoms is given in Eq. (24). 


B. Matrix Elements 

Because the system is homogeneous, the analog of Eq. (16) can be used to express the 
contact operators C on the right side of Eq. (36) in terms of contact density operators at 
the same position r. The wave-vector delta function in Eq. (16) reduces to 6^{K — fco), and 
it can be used to integrate over K. The frequency delta function in Eq. (8) reduces to 

27r(5((EAA — ED)/h ± ca) = 27rd(^u — h/ma^ — fik"^/rn). (37) 

± 

In the sum over ica, only the —u term contributes. 

The expression for the transition rate has been reduced to matrix elements of the contact 
density operator of the form (IT, k\C{r)\kD)- The matrix element is the complex conjugate 
of Eq. (28). The Gaussian integral over K can be evaluated analytically. The sum over hnal 
states of the matrix element in Eq. (36) reduces to 

7 °° k‘^n 

E I I = 64^iVD I dk , (38) 

Before integrating over k, this must be multiplied by the frequency delta function in Eq. (37). 


C. Harmonic Dissociation Rate 


Our hnal result for the harmonic disssociation rate r[^^(a;) in the homogeneous gas can 
be obtained from Eq. (8) by replacing |(/|C'N)P by the right side of Eq. (38), replacing 
the sum of frequency delta functions by the right side of Eq. (37), and then using the 
delta function to integrate over k. The local density approximation can be implemented by 
making the substitution for Nd in Eq. (19). The threshold angular frequency for dissociation 
is h/ma^: the absorption of smaller energy is not enough to break up the dimer. For 
u > h/ma^, the harmonic dissociation rate is 



f bA 

m^al^a [{A + Bo - Sy 


{muj/h — 1/a^)^/^ 
u 


(39) 
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(If the universal dimers are composed of identical bosons, the harmonic dissociation rate is 
given by this same expression.) The harmonic dissociation rate in Eq. (39) has a maximum 
at a; = 2(h/ma^), which is twice the threshold angular frequency. For large frequency, the 
rate decreases very slowly as The dissociation rate is independent of the temperature 

T. This may be surprising at first, but it is related to the fact that the contact of a thermal 
gas of dimers is independent of T. 


D. First-Order Subharmonic Dissociation Rate 


According to Eq. (11), the contribution r 2 ^^(a;) to the subharmonic dissociation rate from 
the hrst-order perturbation in H 2 can be expressed in terms of the harmonic dissociation 
rate in Eq. (39) at twice the frequency. The threshold angular frequency for subharmonic 
dissociation is \{h/ma^). For oj > ^{h/ma^), the transition rate is 







(A + 5o - 



{2mu/n - l/g^)^/^ 

u 


This contribution to the subharmonic dissociation rate has a maximum at a; = %/mb?', which 
is twice the threshold angular frequency. The height of the peak is smaller than that of the 
harmonic dissociation rate by a factor of [6/(A + i?o — .B)]^/4. 


V. SUBHARMONIC ASSOCIATION RATE 

In this section, we calculate the subharmonic association rate in a thermal gas of atoms 
from the second-order perturbation in Hi. We will hnd that this contribution is much larger 
than that from the first-order perturbation in H 2 if B is near a Feshbach resonance. 


A. Initial, Final, and Intermediate States 

We hrst consider a homogeneous gas consisting of Ni atoms of spin state 1 and W atoms 
of spin state 2 in a volume V. If the gas is in thermal equilibrium, the harmonic transition 
rate is given by Eq. (13) with the substitution in Eq. (14), where p = pgas is the density 
matrix for the thermal gas of atoms. In the low-density limit where 3-body and higher- 
body correlations can be neglected, the density matrix pgas can be expressed in terms of the 
density matrix ppair for a pair of atoms in thermal equilibrium, as in Eq. (20). That density 
matrix Ppair is given in Eq. (21). The sum over hnal states reduces to an integral over the 
wave vector of the dimer, as in Eq. (25). 

Once the matrix elements have been reduced to matrix elements in the two-atom sector, 
the sum over intermediate states in Eq. (13) reduces to a sum over atom-pair states and 
dimer states. If the initial state is an atom pair with total energy Eaa, the sum over states 

is 

\^){M ^ [ [ \K',k'){k',K'\ r \k'^){k'^\ 

^mi U JJy {E'p^p^ — Exa)/% ± U} J{E^ — Ep^x)/% U 

where E'^p^ and are given by Eqs. (24) and (26) with primes on the wavenumber variables. 
In the transition rate given by inserting Eq. (17) into Eq. (13), there are four possibilities 


E 
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for the intermediate states \m) and \m') in the amplitude and its complex conjugate: each 
one can be either an atom pair or a dimer. The transition rate can be expressed accordingly 
as the sum of four terms: 


— Taa.AA + rD,D + rAA,D + Td.AA- 
We will calculate each of these terms individually. 


(42) 


B. Matrix Elements 


Because the system is homogeneous, the matrix elements of the contact C in Eq. (13) 
can be expressed in terms of matrix elements of the contact density operator using Eq. (17). 
In addition to the matrix element in Eq. (28) and its complex conjugate, we also need 
the matrix elements of C{r) between atom-pair states and between dimer states. They are 
calculated in Appendix B and given in Eqs. (Bll) and Eqs. (B13): 


{K',k'\C{r)\K,k) 

{Wr)\k^) 


IbTT^a^ 

(1 -|- m/c')(l — iak) ’ 
8n/a. 


(43a) 

(43b) 


The integrals over the total wave vectors of the intermediate states and over the wave 
vector of the hnal-state dimer can be evaluated using the delta functions in Eq. (17). The 
Gaussian integral over the total wave vector of the initial atom-pair state can then be 
evaluated analytically. The frequency delta function reduces to 

27r(5((i7D — Eaa)/^ ± 2a;) = 2ttS(2u — h/ma^ — /m). (44) 

± 


In the sum over ±2a;, only the -|-2a; term contributes. 


1. Intermediate atom-pair states 


The contribution from intermediate atom-pair states to the factor in the transition rate 
involving matrix elements reduces to 


EE 

i / 


E 


(/|C'|m)(m|C'|i) 


A^lA^2(i|Ppair|i) 


8192y27r 


N 1 N 2 

~l?d V~ 


/““, A;^exp(— 


r*cxD 

dk 


k 


Uo {k"^ + 1/a?){k'‘^ — k‘^ + muj/h)\ 


<45) 


Before integrating over X, this must be multiplied by the frequency delta function in Eq. (44). 

The threshold angular frequency for subharmonic association is \{?i/md?). For u > 
h/md?', the integral over k' in Eq. (45) has a pole on the integration contour. In this region 
of a;, this contribution to the transition rate is a subleading correction of order b'^ to the 
harmonic transition rate of order in Eq. (30). We therefore consider only the frequency 
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interval \{h/ma^) < u < h/ma^, where this contribution to the transition rate is leading 
order in b. In this region of ca, the integral over k' in Eq. (45) is 



k'^ 

(A;'2 + l/a2)(A;'2_ 


k? + moj/h) 


Tia 

2(1 + a^muj/h — k'^) 


(46) 


The frequency delta function in Eq. (44) can be used to evaluate the integral over k in 
Eq. (45). The resulting contribution to the transition rate is 


rAA,AA = 




hA 


(A + Bo-B)^ 

where k,{2u) is given in Eq. (33). 


N 1 N 2 k{2uj) exp ( — j3h^K^{2uj)/m) 

^ a;(l + a/1 — mua?/hY 


(47) 


2. Intermediate dimer states 


The contribution from intermediate dimer states to the factor in the transition rate 
involving matrix elements reduces to 


* / 


E 


{f\C\m){m\C\i) 




NiN2{i\p^irLr\i) 


[A . 


V Jq {k"^ + 1/a^){k‘^ + 1/a? — moj/hy 


(48) 


Before integrating over k, this must be multiplied by the frequency delta function in Eq. (44). 
The frequency delta function can be used to evaluate the integral over k in Eq. (48). The 
resulting contribution to the transition rate is 


V2Ty\% ( bA 

\{A + Bo-By 


4 


iViiV2 k,(2u) 


exp (-/3h\^(2a;)/m) . 


(49) 


3. Interference between Atom-Pair and Dimer States 


The contribution to the factor in the transition rate involving matrix elements from 
intermediate atom-pair states in the amplitude and from intermediate dimer states in its 
complex conjugate reduces to 


LEE „ . +,, E —— 


f rn 




-8192V2x^ ’"dj 


r*cxD 

dk- 


k'^ exp{—/3fyk'^/m 


V Jq (A;2l/a2)(A;2l/a2 - ca/h) 

f°° t k'"^ 

Jo {k'‘^1/dy {k"^ — k'^mu/h) ^ ^ 
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Before integrating over k, this must be multiplied by the frequency delta function in Eq. (44). 
The integral over k is given in Eq. (46). The frequency delta function can be used to evaluate 
the integral over k in Eq. (48). The resulting contribution to the transition rate is 

^ ^ 6A iViiV 2 fi:(2a;) exp (—/dh^K^(2a;)/m) 

\(A + 5o — 5)^/ V — mua'^/h) 

The contribution rD,AA from intermediate dimer states in the amplitude and intermediate 
atom-pair states in its complex conjugate is the same as Eq. (51) for frequencies in the range 
\{h/ma^) < Lo < h/ma^. The contributions Taa.d = rD,AA are negative, because there is 
destructive interference between atom-pair and dimer intermediate states. 


C. Total subharmonic transition rate 


The total subharmonic transition rate in Eq. (42) from the second-order perturbation in 
Hi is given by adding Eqs. (47) and (49) and twice Eq. (51). The local density approximation 
can be implemented by making the substitution for N 1 N 2 /V in Eq. (18). For frequencies in 
the range \{h/ma^) < u < h/ma^, the subharmonic transition rate is 


p(iii) 

2 


(ca) 


\/2h^a / 6A 
4m2a4g l^(A + 5o-5)2 



X exp {-(3h\\2uj)/m) 

u ^ ' 


1 -|- a /1 — mua^/h 


muja?/h j 


.(52) 


where k{2uj) is given by Eq. (33). (The corresponding result for identical bosons can be 
obtained by replacing ni{r)n 2 {r) by n^(r)/2, where n{r) is the local number density.) 

If ksT k?/ma^^ the subharmonic association rate in Eq. (52) has a narrow peak with a 
maximum when uj is above the threshold ^{h/ma^) by approximately ksT/Ah. In the region 
near the threshold and the peak, the largest contribution comes from intermediate dimer 
states. The intermediate atom-pair states give a contribution that is smaller at threshold 
by a factor of (3 — 2\/2)/S ~ 0.021. The cross terms give a negative contribution that is 
smaller at threshold by a factor of (2 — a/2)/2 ^ 0.29. The subharmonic association rate in 
Eq. (52) is much smaller than the harmonic association rate in Eq. (30). The ratio of their 
maximum values is 


T 

T 


( 1 , 1 ) 

2,max 

l,max 


2.914 


6A 


(A + 5o - 





1.21 


ksTma 



(53) 


(2) 

The contribution T^ to the subharmonic transition rate from the hrst-order perturbation 
in H 2 is given in Eq. (32). Near the subharmonic threshold frequency, differs from T^^^ 
by a factor of 11.6 (a/abg)^[A/(A + Bq — If B is near the Feshbach resonance, T^^^ 

is much smaller. It is therefore unnecessary to consider interference between the hrst-order 
perturbation in Hi and the second-order perturbation in H 2 - 
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VI. SUBHARMONIC DISSOCIATION RATE 


In this section, we calculate the subharmonic dissociation rate in a thermal gas of dimers 
from the second-order perturbation in Hi. We hrst consider a homogeneous gas of N-q dimers 
in a volume V in thermal equilibrium at temperature T. The subharmonic transition rate is 
given by Eq. (13) with Eq. (17) inserted and with |i)(i| replaced by the density matrix pgas 
for the thermal gas of dimers. In the low-density limit where correlations between dimers 
can be neglected, pgas can be expressed in terms of the density matrix pdimer for a dimer in 
thermal equilibrium, as in Eq. (34). The density matrix pdimer is given in Eq. (35). The sum 
over hnal states reduces to an integral over center-of-mass wave vector K and relative wave 
vector /c of a pair of atoms, as in Eq. (36). The frequency delta function reduces to 

2'K5(y{ExA — Ejy)/h ± 2oj) = 2 tt6(2u — h/ma'^ — /m). (54) 

± 


In the sum over ±2ca, only the —2uj term contributes. In the sums over intermediate states 
in the amplitude and its complex conjugate, both intermediate states can be either an atom 
pair or a dimer. The calculation of the individual contributions proceeds in the same way 
as for the association rate. They can be obtained from those in Eqs. (47), (49), and (51) by 
replacing N 1 N 2 /V by Ad and replacing A^exp (^—l3h‘^K^{2u)/m) by v^/4. 

In the local density approximation, the factor Ad is replaced by J d^rn-£,{r). For fre¬ 
quencies in the range \{h/ma^) < u < h/ma^, the subharmonic dissociation rate is 


■p(iii) 

J- 2 


(cn) 


h^a 




( 1 

^ \^ + a/1 — mud?/h 



(55) 


where k{2u) is given by Eq. (33). (If the universal dimers are composed of identical bosons, 
the subharmonic dissociation rate is given by this same expression.) Like the harmonic 
dissociation rate in Eq. (39), the subharmonic dissociation rate in Eq. (55) is independent 
of the temperature T. The subharmonic dissociation rate has a maximum at an angular 
frequency u that is above the threshold \{h/md?) by approximately 0.082(h/ma^). The 
subharmonic dissociation rate in Eq. (52) is much smaller than the harmonic dissociation 
rate in Eq. (39). The ratio of their maximum values is 


T 

T 


(1,1) 

2, max 

l,max 


1.39 


5A 


(A + Ao - Bf 



(56) 


i<y\ 

The contribution T^ to the subharmonic transition rate from the first-order perturbation 
in H 2 is given in Eq. (40). Near the subharmonic threshold frequency, differs from T^^^ 
by a factor of 11.6 (a/abg)^[A/(A + Bq — B)Y. If B is near the Feshbach resonance, T^^^ 
is much smaller. It is therefore unnecessary to consider interference between the hrst-order 
perturbation in Hi and the second-order perturbation in H 2 . 
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VII. APPLICATION TO ^LI ATOMS 


An experiment on the association of ^Li atoms into universal dimers using a modulated 
magnetic field was performed at Rice University by Dyke, Pollack, and Hulet [8]. The 
^Li atoms were in \F = l^mp = 1) state. The scattering length was controlled using 
the Feshbach resonance at Bq = 737.7 G. The other parameters in the expression for the 
scattering length in Eq. (2), as determined in Ref. [8], are Obg = —20.0 oq and A = —174 G. 
The bias magnetic field was set to B = 734.5 G. The measurements were carried out for 
both a Bose-Einstein condensate and a thermal gas of ^Li atoms. Magneto-association into 
dimers was observed through the loss of atoms, presumably from inelastic collisions of the 
dimers with atoms. 

In the experiment on the BEG of ’^Li atoms in Ref. [8], they observed a narrow loss 
resonance as a function of the frequency, with the fraction of atoms remaining decreasing 
almost to 0. The position of this resonance was used to measure the binding energy of the 
dimer to be h(450 kHz), corresponding to a ~ 1100 Oq. 

The experiment on the thermal gas of ^Li atoms in Ref. [8] was carried out at three 
combinations of the amplitude b of the modulated magnetic held and the temperature T: 
(b, T) = (0.57 G, 3 /iK), (0.14 G, 3 /iK), and (0.57 G, 10 fj,K). The duration of modulation was 
in the range from 25 fis to 500 /is, but its value was not specihed for each individual set (6, T). 
The fraction of atoms remaining after the modulation time was measured as a function of 
frequency of the oscillating held. The data were ht to convolutions of Lorentzians with the 
thermal Boltzmann distributions. For each set {b,T), there was a harmonic peak just above 
ud- a subharmonic peak just above was evident only for (0.57G, 3/iK). For these 

values of {b,T), the minimum fractions remaining were about 0.3 in the harmonic peak and 
about 0.5 in the subharmonic peak. 

Since the number of atoms in the thermal clouds and the modulation times for each set 
of {b,T) were not specihed in Ref. [8], we are unable to make quantitative comparisons with 
our theoretical results. In Fig. 1, we show our results for the association rates as functions 
of the angular frequency u for the three sets of values of {b,T) for which the atom loss was 
measured in Ref. [8]. The ratio of ksT/h to the binding frequency ud = fi/ma? of the dimer 
is 0.15 and 0.49 at the temperatures 3/iK and 10/iK, respectively. Thus the condition 
ksT h^/ma^ is much better satished at 3/iK. In Fig. 1, the curves for uj > ud are 
the harmonic association rates F^^^(a;) in Eq. (30) with J d^rnin 2 replaced by J d^rin? /2. 
The curves for oj < uo are the subharmonic association rates F 2 ^’^^(a;) in Eq. (52) with 
J d?rnin 2 replaced by j d^rri^ /2. The subharmonic association rates F 2 ^^(a;) are completely 
negligible: their peak values are smaller than those for by factors of about 3 x 10“^. 
The association rates in Fig. 1 are divided by j d^rn^{r) = N (n) to obtain rates F/iV(n) 
that do not depend on the number of atoms. The angular frequency is normalized to the 
angular binding frequency uo = h/md?' of the dimer. The heights of the harmonic and 
subharmonic peaks for (5, T) = (0.57 G, 3 /xK) are larger than those for (0.14 G, 3 /iK) by the 
ratio of the values of 5^, which is 12.8. The maxima of the harmonic association rates are 
at angular frequencies u that are above the threshold wd by approximately ksT/2h. The 
maxima of the subharmonic association rates are at angular frequencies that are above the 
threshold by approximately ksT/Ah. The ratios of the maximum of the harmonic peak 
to the maximum of the subharmonic peak are 11.4, 190, and 13.8 for (&, T) = (0.57 G,3/iK), 
(0.14 G, 3/iK), and (0.57G, 10/iK), respectively. They can be compared to the ratios 11.9, 
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(O/CO^ 


FIG. 1: Association rate r/A^(re) for a thermal gas of ^Li atoms as a function of the modulation 
frequency of the magnetic field. The angular frequency a; is normalized to the angular binding 
frequency of the dimer. The bias field is .B = 734.5 G, which corresponds to a ~ 1100 uq. The 
three sets of curves are for different values of the modulation amplitude h and the temperature T: 
(6, T)=(0.57 G, 3 ^K) (highest peak value, red), (0.57 G, 10/xK) (green), and (0.14 G, 3/rK) (lowest 
peak value, blue). The association rates in the region ui < ujd have been multiplied by 10 to make 
the subharmonic transitions more visible. 

198, and 24.1 predicted using Eq. (53). 

Our results for dissociation rates in a thermal gas of dimers as functions of the angular 
frequency uj are illustrated in Fig. 1 using the two values of b for which the atom loss 
was measured in Ref. [8]. The dissociation rates in Fig. 1 are divided by the number N^, 
of dimers to obtain rates T/No that do not depend on No- The angular frequency u is 
normalized to the binding angular frequency ud oi the dimer. The curves for u > ud 
are the harmonic dissociation rates r^^^(a;) in Eq. (39). The curves for u < ojd are the 
subharmonic dissociation rates r 2 ^’^^(a;) in Eq. (55). The subharmonic association rates 
r 2 ^^(a;) are completely negligible: their peak values are smaller than those for by 

factors of about 3 x 10“®. The heights of the harmonic and subharmonic peaks for b = 0.57 G 
are larger than those for b = 0.14 G by the ratio of the values of 6^, which is 12.8. The 
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FIG. 2: Dissociation rate T/Nd for a thermal gas of ^Li dimers as a function of the angular 
frequency oj of the modulated magnetic field. The angular frequency u is normalized to the 
binding angular frequency ujd of the dimer. The bias field is .B = 734.5 G, which corresponds to 
a ~ 1100 ag. The dissociation rate is independent of the temperature. The two sets of curves are 
for different values of the modulation amplitude b: 0.57 G (higher peak value, red) and 0.14 G 
(lower peak value, blue). The dissociation rates in the region oj < ujd have been multiplied by 10 
to make the subharmonic transitions more visible. 

maxima of the harmonic association rates are at twice the threshold angular frequency 
Ti/ma^. For large cu, the harmonic rates decrease very slowly as The maxima of the 

subharmonic dissociation rates are above the threshold angular frequency \{h/ma?‘) by only 
about 0.08 h/ma?. The ratio of the maximum of the harmonic peak to the maximum of 
the subharmonic peak is given in Eq. (56). The ratios for b = 0.57 G and b = 0.14 G are 
approximately 21 and 340, respectively. 


VIII. COMPARISONS AND SUMMARY 

In this section, we describe previous theoretical treatments of the association of atoms 
into dimers using a modulated magnetic held. We then summarize our results on association 


19 









and dissociation and compare the association results with those from the other approaches. 

A theoretical treatment of the association of atoms into dimers using an oscillating mag¬ 
netic held was presented by Hanna, Kohler, and Burnett in 2007 [10]. They used a two- 
channel model for the two-atom system, with one channel consisting of atom pairs interacting 
through a short-range separable potential and a second channel consisting of a single discrete 
molecular state. The dimer is an eigenstate of the coupled-channel problem. The modula¬ 
tion of the magnetic held was taken into account through the sinusoidal oscillation of the 
energy of the discrete molecular state. They solved the time-dependent Schrodinger equa¬ 
tion for the two-channel model numerically to obtain the probability for association into the 
dimer as a function of time. They studied the dependence of the association probability for 
a homogeneous thermal gas on the frequency and amplitude of the oscillating held and on 
the temperature and density of the gas. Association into dimers in a Bose-Einstein conden¬ 
sate was treated in a completely diherent way by solving numerically an integro-diherential 
equation for the mean-held of a homogeneous BEC. The results for association probabili¬ 
ties in Ref. [10] are all completely numerical. They can be applied quantitatitively only to 
homogeneous systems of ®^Rb and ^^^Cs atoms under the specihc conditions considered in 
the paper. In order to use their methods to predict the association probability for any other 
conditions, such as diherent atoms, other oscillation parameters, diherent temperature or 
density for the homogeneous gas, or a trapped gas with variable density, it would be neces¬ 
sary to solve their equations numerically for each set of conditions. Subharmonic transitions 
were not considered in Ref. [10]. 

Brouard and Plata have recently presented a diherent theoretical treatment of the as¬ 
sociation of atoms into dimers using an oscillating magnetic held [11]. They also used a 
two-channel model for the two-atom system, with one channel consisting of a continuum 
of positive-energy atom pair states and the second channel consisting of a single discrete 
molecular state. The dimer is an eigenstate of the coupled-channel problem. As in Ref. [10], 
the modulation of the magnetic held was taken into account through the sinusoidal oscilla¬ 
tion of the energy of the discrete molecular state. However in Ref. [11], a time-dependent 
unitary transformation was used to move the time-dependence into oh-diagonal terms be¬ 
tween the dimer state and the atom-pair states. They showed that for frequencies near 
the harmonic resonance, the dynamics in the transformed frame can be approximated by a 
time-independent Hamiltonian whose entries are given analytically as functions of the oscil¬ 
lation parameters and the Feshbach resonance parameters. Subharmonic transitions in an 
appropriate transformed frame are described by a diherent time-independent Hamiltonian 
whose entries are given analytically. Some qualitative aspects of the association process were 
deduced from these ehective Hamiltonians. However the results in Ref. [11] for association 
probabilities in thermal gases of ®^Rb atoms and in Bose-Einstein condensates of ®®Rb atoms 
are completely numerical. 

In this paper, we have applied the new approach to this problem that was introduced in 
Ref. [12]. It was based on the realization that the leading ehect of an oscillating magnetic held 
near a Feshbach resonance can be treated as a time-dependent perturbation proportional 
to the contact operator C. In Appendix A, we presented a quantum held theory argument 
that the perturbation proportional to C can also be used beyond hrst order. Fermi’s Golden 
Rule is used to obtain general expressions for transition rates in terms of transition matrix 
elements of C. Our general formula for the harmonic transition rate in Eq. (8) comes 

from the hrst-order perturbation in C and was obtained previously in Ref. [12]. There is a 
contribution r 2 ^^(a;) to the subharmonic transition rate from the hrst-order perturbation in 
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C that is related in a simple way to the harmonic rate and is given in Eq. (11). However 
the second-order perturbation in C gives another contribution r 2 ^’^^(a;) to the subharmonic 
transition rate that is given in Eq. (13). Near the Feshbach resonance, is larger than 
by a factor of (a/obg)^. 

For a homogeneous system, our general expressions for the transition rates can be sim- 
plihed by expressing the transition matrix elements of the contact operator in terms of 
transition matrix elements of the contact density operator C. The harmonic transition rate 
in Eq. (8) can be simplihed by inserting Eq. (16). The subharmonic transition rate 
in Eq. (13) can be simplihed by inserting Eq. (17). For a nonhomogeneous system 
in the local density approximation, these simplihcations can hrst be used to calculate the 
transition rates for the homogeneous system. Substitutions such as those in Eqs. (18) and 
(19) can then be used to obtain the transitions rate for the nonhomogenous system. 

To obtain association rates in a thermal gas of atoms and dissociation rates in a thermal 
gas of dimers, we hrst exploited the low density to reduce the transition matrix elements of C 
in the thermal gas to transition matrix elements of C in the two-body problem. Those matrix 
elements were calculated in Appendix B using the quantum held theory formulation of the 
problem of atoms with zero-range interactions. In the two-atom sector, this is equivalent 
to a single-channel model for atoms with large scattering length, with the dimer arising 
dynamically as a bound state. This allowed us to calculate the matrix elements of the 
contact density analytically. 

Our hnal results for the harmonic and subharmonic association rates in a thermal gas of 
atoms are given in Eqs. (30) and (52). Our hnal results for the harmonic and subharmonic 
disssociation rates in a thermal gas of dimers are given in Eqs. (39) and (55). These results 
are analytic functions of all the relevant parameters: the oscillation parameters ca, b, and 
a or B, the Feshbach resonance parameters Obg, Bq, and A, and the temperature T. The 
association rates in a thermal gas of fermions with two spin states depend on the local 
number densities ni{r) and n 2 {r) only through the multiplicative factor J d^rnin 2 - The 
dissociation rates in a thermal gas of dimers depend on the local number density n-D{r) 
only through the multiplicative factor f d^rriB- Our analytic results should be useful for 
analyzing experiments on association into and dissociation of dimers. They should also be 
useful for designing experiments that optimize the number of dimers created or destroyed 
by the modulated magnetic held. For a thermal gas of atoms with ksT b?/md?^ the 
maximum in the harmonic association rate is at an angular frequency u that is above the 
threshold Ti/ma^ by approximately kBT/2h. The maximum in the subharmonic association 
rate is at an angular frequency that is above the threshold |(h/ma^) by approximately 
ksT/Ah. For a thermal gas of dimers, the maximum in the harmonic dissociation rate is at 
an angular frequency u that is approximately twice the threshold Ti/ma^. The maximum 
in the subharmonic dissociation rate is at an angular frequency that is above the threshold 
\{h/md?) by approximately 0.08 {h/md?). 

Our general results for the harmonic and subharmonic transition rates in terms of matrix 
elements of the contact density operator can also be applied to other systems. An analytic 
result for the association rate in a dilute Bose-Einstein condensate of identical bosons was 
given in Ref. [12]. It should also be possible to obtain analytic results for superfluids of 
fermions with two spin states at zero temperature, including the dissociation rate of dimers 
in the BEC limit and the dissociation rate of Cooper pairs in the BCS limit. The dissociation 
rate of paired fermions in the unitary limit is more challenging, but it is an important problem 
because it would allow the first direct measurements of the gap for the unitary Fermi gas. 
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Appendix A: Quantum Field Theory Derivation of Perturbing Hamiltonian 


Particles with a scattering length a that is large compared to the range tq of their in¬ 
teractions can be described by a local quantum held theory. For a fermion with two spin 
states, there are two fermionic quantum helds 'ipi and ip 2 - The interactions of the quantum 
held theory are made local by taking the zero-range limit at the expense of introducing 
an ultraviolet cutoh A on the momenta of virtual particles. The interaction Hamiltonian 
density is 

■Hint = (Al) 

where Aq is the bare coupling constant. If h is set to 1, Aq has dimensions of length. The 
held theory describes particles with scattering length a if the bare coupling constant is 


Ao — 


dTT 

1/a — 2A/7r 


(A2) 


Matrix elements of the operator diverge as A^ as the cutoh is increased to cx). 

Since Aq scales as 1/A, matrix elements of the interaction Hamiltonian density in Eq. (Al) 
therefore diverge as A. In matrix elements of the complete Hamiltonian density, the diver¬ 
gence is cancelled by a corresponding divergence in matrix elements of the kinetic energy 
density. In matrix elements of the operator subleading terms that diverge as A 

give hnite contributions to the energy density. The contact density operator in the quantum 
held theory is [15] 

C = (A3) 

This operator has hnite matrix elements, because the divergence in the matrix element of 
'ip\'ipl'ip 2 'ipi proportional to A^ is compensated by a factor of 1/A^ from Ag. 

The local quantum held theory can describe particles with a time-dependent scattering 
length a(t) provided the time scale a/a is large compared to the time scale mrl/U set by the 
range. In the interaction Hamiltonian density in Eq. (Al), the time-dependent bare coupling 
constant Ao(t) is obtained by replacing a in Eq. (A2) by a(t). If the time dependence consists 
of small deviations in the inverse scattering length from some value 1/a, the bare coupling 
constant can be expanded around the corresponding value Aq: 


Ao(^) 


dvr \a{t) 


(— 

aj (47r)2 \o(^) 



(A4) 


When this is inserted into the interaction Hamiltonian density in Eq. (Al), the term linear 
in l/a{t) is proportional to the contact density operator in Eq. (A3). The term quadratic 
in l/a{t) is suppressed by 1/A from the additional power of Ag. The higher order terms are 
even more highly suppressed. Thus the interaction Hamiltonian density in the zero-range 
limit can be reduced to 


'Hint(^) 


m 


(— 

4:71 m \a(t) 



(A5) 
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Appendix B: Matrix Elements of the Contact density operator 


The field theoretic definition of the contact density operator in Eq. (A3) can be expressed 
as 

C{r) = 0t(r)0(r), (Bl) 

where the contact field is a local operator that annihilates two atoms at a point. 

In the case a > 0, 4>{r) has a nonzero amplitude to annihilate a dimer, so it can also be 
referred to as the dimer field. The transition matrix element of the contact density operator 
can be expressed as 

{f\C{r)\i) = ^(/|0i(r)|n)(n|0(r)|i). (B2) 

n 

A complete set of states = 1 has been inserted between and fi. If only one term 

in the sum is nonzero, the matrix element factors into a matrix element of that involves 
the initial state and a matrix element of that involves the final state. 

In order to calculate magneto-transition rates in a thermal gas of atoms or dimers, one 
needs to calculate transition matrix elements of the contact density operator C{r) between 
two-atom states, which are either a pair of unbound atoms or a dimer. We will calculate these 
matrix elements in the Zero-range Model defined by the interaction Hamiltonian density in 
Eq. (Al). The Feynman rules for the atom propagator and the 2-atom-to-2-atom vertex are 
specified in the appendix of Ref. [16]. The 2-atom-to-molecule coupling constant go should 
be set to 0. Using these Feynman rules, the calculation of transition matrix elements of the 
contact density operator can be reduced to evaluating Feynman diagrams. 

The transition amplitude A{Ecm) is the amplitude for the transition between a pair of 
atoms in the asymptotic past and a pair of atoms in the asymptotic future. It is a function 
only of the energy Eon of the pair of atoms in their center-of-mass frame: 

Ecra = E — K'^/Am, (B3) 

where E is their total energy and K is their total momentum. (We set h = 1 in this 
Appendix.) The transition amplitude can be calculated by solving the Lippmann-Schwinger 
equation shown in Figure 3: 


iA{Eca,) = -i{\o/m) + {\o/m)I{Eca,)A{Ea„)- 


(B4) 


The loop integral I{Ecm) in the last diagram in Figure 3 is 


I{Ecm) 


II f 1 

fn J.J - mEcrn - ie ’ 


(B5) 


where q is the loop momentum. Using the expression for the bare coupling constant in 
Eq. (A2), the solution can be expressed as 


A(F'cm) 


47r/m 

— 1/a -|- \/—mEcra — A 


(B6) 


This amplitude has a pole in the energy ai E = K‘^/4m — 1/ma?. The residue of the pole 
is —Zy), where 

Zd = Sn/m'^a. (B7) 
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FIG. 3: The Lippmann-Schwinger equation for the transition amplitude A(Ecui)- The blob repre¬ 
sents iA{Ecra)- The vertex factor is —iXojm. 

The standard Feynman rnles can be nsed to calculate matrix elements of local operators 
between states in the asymptotic past and states in the asymptotic future. However we need 
matrix elements between initial and hnal states at the same time. We will use the Feynman 
rules to calculate matrix elements of the contact held operator between the vacuum and 
two-atom states in the asymptotic future. We will then calculate matrix elements of the 
contact density operator between two-atom states in the asymptotic future by expressing 
them in terms of matrix elements of <p^ and 0. 


1. Vacuum-to-pair matrix element of 



F 



FIG. 4: Feynman diagrams for the matrix element of 0^(r) between the vacuum and the atom-pair 
state \K,k). The open dot represents the operator, whose Feynman rule is Aq. The residue of 
the pole in the energy in the second diagram can be used to determine the matrix element of 0l (r) 
between the vacuum and the dimer state {ku). 

The matrix element of between the vacuum |0) and an atom-pair state \K, k) with 
total momentum K and relative momentum k can be calculated from the Feynman diagrams 
in Fig. 4. The atoms are on their energy shells with total energy E = /Am + k‘^/m. The 
matrix element is 

(K,fc|0t(r)|O) = + 

= -mH(Ecm). (B8) 

The energy E^m of the atom pair in their center-of-mass frame depends only on their relative 
momentum: 

Ecra = k^/m. (B9) 


2. Vacuum-to-dimer matrix element of 

The matrix element of 4>^{r) between the vacuum |0) and a dimer state I^d) with momen¬ 
tum fen can be calculated from second Feynman diagram in Fig. 4. That diagram has a pole 
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in the total energy E of the hnal-state atoms, which must be off their energy shells. At the 
pole, the center-of-mass energy is equal to the binding energy of the dimer: = —1/mo?. 

The residue of the pole is the product of the desired matrix element and where Zd is 
the residue factor given in Eq. (B7). The matrix element is therefore 

(fcD|0^(r)|O) = ^Ao/(Ee^)(-ZD)Z-'/^ 

= a/Stt/o. (bio) 


3. Pair-to-pair matrix element of 

The matrix element of 010 between the atom-pair states \K,k) and \K\k') in the 
asymptotic future can be calculated by inserting a complete set of states between 0l and 0, 
as in Eq. (B2) . Since the operator 0 annihilates the initial-state atoms, the only term that 
contributes is the vacuum state. The matrix element factors into the vacuum-to-atom-pair 
matrix element and its complex conjugate: 

{K\k’\4,Hr)4,{r)\K,k) = (if', fc'1^(r)|0)(0|y (r)|iC, fe>, 

= (Bll) 

where Ecm is given in Eq. (B3) and E'^^ is the same expression with k replace by k'. 

The expression for this matrix element given in Ref. [16] is incorrect: the factor A*(Ecm) 
was not complex conjugated. It is easy to see that this is incorrect by setting the hnal 
state equal to the initial state: K = K, k = k. Since the operator is hermitian, the 
matrix element must be real. This condition is satisfied by Eq. (Bll). The error made in 
Ref. [16] was that the matrix element was calculated not between atom-pair states at the 
same time, but between an atom-pair state in the asymptotic past and an atom-pair state 
in the asymptotic future. 


4. Pair-to-dimer matrix element of 0^0 

The matrix element of 0l0 between the atom-pair state \K, k) in the asymptotic future 
and the dimer state Ifc^) in the asymptotic future can be calculated by inserting a complete 
set of states between 0f and 0, as in Eq. (B2) . Since the operator 0 annihilates the initial- 
state atoms, the only term that contributes is the vacuum state. The matrix element factors 
into a vacuum-to-dimer matrix element and the complex conjugate of a vacuum-to-atom- 
pair matrix element: 

(fcD|0Hr)0(r)|R:,fc) = (fcD|0Hr)|O)(O|0(r)|R:,fc), 

= -m^/8n/aA*{Ecm), (B12) 


where Ecm is given by Eq. (B3). 


5. Dimer-to-dimer matrix element of 0^0 

The matrix element of 0^0 between the dimer states I^d) and |fcp) in the asymptotic 
future can be calculated by inserting a complete set of states between 0i and 0, as in Eq. (B2). 
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Since the operator 0 annihilates the initial-state dimer, the only term that contributes is 
the vacuum state. The matrix element factors into a vacuum-to-dimer matrix element and 
its complex conjugate: 


{k'^\(j)\r)(j){r)\kY)) = {k'^\4>\r)\0) {0\(j){r)\kj)), 
= Sn/a. 


(B13) 


[1] S.T. Thompson, E. Hodby, and C.E. Wieman, Phys. Rev. Lett. 95, 190404 (2005) [cond- 
mat/0505567]. 

[2] S.B. Papp and C.E. Wieman, Phys. Rev. Lett. 97, 180404 (2006) [cond-mat/0607667]. 

[3] C. Weber, G. Barontini, J. Catani, G. Thalhammer, M. Inguscio, and F. Minardi, Phys. Rev. 
A 78, 061601 (R) (2008) [arXiv:0808.4077]. 

[4] A.D. Lange, K. Pilch, A. Prantner, F. Ferlaino, B. Engeser, H.-C. Naegerl, R. Grimm, and 
G. Chin, Phys. Rev. A 79, 013622 (2009) [arXiv:0810.5503]. 

[5] N. Gross, Z. Shotan, O. Machtey, S. Kokkelmans, and L. Khaykovich, Comptes Rendus 
Physique 12, 4 (2011) [arXiv:1009.0926]. 

[6] S.E. Pollack, D. Dries, R.G. Hulet, K.M.F. Magalhaes, E.A.L. Henn, E.R.F. Ramos, 
M.A. Caracanhas, and V.S. Bagnato, Phys. Rev. A 82, 020701(R) (2010) [arXiv: 1004.2887]. 

[7] O. Machtey, Z. Shotan, N. Gross, and L. Khaykovich Phys. Rev. Lett. 108, 210406 (2012) 
[arXiv: 1201.2396]. 

[8] P. Dyke, S.E. Pollack, and R.G. Hulet, Phys. Rev. A 88, 023625 (2013) [arXiv:1302.0281]. 

[9] D.H. Smith, arXiv: 1503.02688. 

[10] T.M. Hanna, T. Koehler, and K. Burnett, Phys. Rev. A 75, 013606 (2007) [arXiv:cond- 
mat/0609725]. 

[11] S. Brouard and J. Plata, J. Phys. B: At. Mol. Opt. Phys. 48, 065002 (2015) [arXiv: 1503.01700]. 

[12] G. Langmack, D.H. Smith, and E. Braaten, Phys. Rev. Lett. 114, 103002 (2015) 
[arXiv: 1406.7313]. 

[13] S. Tan, Ann. Phys. 323, 2971 (2008) [cond-mat/0508320]. 

[14] E. Braaten and H. -W. Hammer, J. Phys. B 46, 215203 (2013) [arXiv:1302.5617]. 

[15] E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 (2008) [arXiv:0803.1125]. 

[16] E. Braaten, D. Kang, L. Platter, Phys. Rev. A. 78, 053606 (2008) [arXiv:0806.2277]. 


26 



